SOLITON DYNAMICS FOR A GENERAL 
CLASS OF SCHRODINGER EQUATIONS 



RAFFAELLA SERVADEI* AND MARCO SQUASSINAt 

Abstract. The soliton dynamics for a general class of nonlinear focusing Schrodinger problems 
in presence of non-constant external (local and nonlocal) potentials is studied by taking as initial 
datum the ground state solution of an associated autonomous elliptic equation. 



1. Introduction and main result 

1.1. Introduction. The aim of this paper is to study a general class of scalar and vectorial 
Schrodinger equations in presence of local and nonlocal potentials, modelling an electric and mag- 
netic field and a Newtonian type interaction, respectively. This class of problems includes various 
physically meaningful particular cases, that will be individually described in details later in this 
section. In fact, we would also like to discuss the latest developments available in literature for 
this kind of issue, particularly when approached via the technique initiated by the 2000 work of 
R. Jerrard and J. Bronski [BJ]. More precisely, let m > 1, N > 1, < p < 2/iV, e > and let 

(1.1) V : R N — > R, A : R N -► R N , $ : R^ -» R, 

be C 3 (M N ) functions satisfying suitable assumptions that will be stated in the following. Then, if i 
denotes the complex imaginary unit, consider the Schrodinger equation 

■iedtd + L A (i + V(x)Ci = \(s\f(i + * \Q) CI in R N x (0, oo), 
(S) {(i(x,0)=(i(x) inR N , 

J = 1, • • • ,m, 

where Q = . . . , Q 1 ) : R N x R + — ► C' m is the unknown, the magnetic operator La is defined as 

L A C := -yAC - U(x) ■ VC + \\A(x)\ 2 C - ^div x A(x)(, 
the convolution is denoted by ($ * v)(x) := J <&(x — y)v(y)dy, and 

m m 

icif :=a,-ic j f p +X>dct + Vr\ ICg-^ICf+IJwylC*! 2 , 
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for some nonnegative constants ai, fa, 7^ , such that 7^ = 7^ and uii 



for all i,j 



, m. 



By rescaling problem (S) with 4> £ (x,t) = ( £ (ex,et), we reach the following system, where e appears 
now only in the arguments of the potentials V, A and $ 



(P) 



with 



(1.2) 



-id t 4 + L A 4 + V{ex)fr 

4(x,o) = 4(x) 

lj = 1,... ,m, 



m 



m 



a A? 



X (0, 00) 



6 m) 



xR + ^C m and 
~A(f) - jA{ex) ■ V(f> + -\A(ex) 



— div x A(ex) 



* are an electric and magnetic 
and can be thought (and identified) 



As we have already recalled, here V : 1^ — > K and ^4 : 
potentials respectively. The magnetic field £>is-B = Vx^4inI 
in general dimension as a 2-form H s of coefficients (diAj — djAi). We will keep using the notation 
B = V x A in any dimension iV. 

We point out that the general Schrodinger problem (5*) we aim to investigate contains, as particular 
cases, the following physically meaningful situations. 



Class I. If m = l, A = 0, fa 



jij = and ctj = 1, one finds: 



'ied t Ce + 4 AC e - V(x)Ce + \Ce\ 2p Ce = in R N X (0, oo), 
Xe(x,0) = Co(x) in 



This is the classical Schrodinger equation with a spatial potential. For general results about local and 
global existence of solutions, regularity, orbital stability and instability, we refer the reader to [Ca] 
and to the references therein. From the point of view of the semi-classical analysis of standing 
wave solutions £ £ (x,t) = u £ (x)e~ lEt for E € M, the Schrodinger equation reduces to a semi-linear 
elliptic equation. In the last few years a huge literature has developed starting from the celebrated 
paper by Floer and Weinstein [FW] (see the monograph [AM] by Ambrosetti and Malchiodi and 
references therein). Concerning the soliton (or, equivalently, point-particle) dynamics, that is the 
study of the qualitative behaviour of the solutions of this equation by choosing as initial datum 
a suitably rescaled ground state solution of an associated elliptic problem, we refer e.g. to the 
works [BJ, FGJS, GSS, Ke] and to the recent monograph [Car] (see also e.g. [KN, KM] for works 
in the mathematical physics community). Very recently, in [BGM], Benci, Ghimenti and Micheletti 
provided the first result on the soliton dynamics with uniform global estimates in time. 



Class II. If m = 1, Pj 



UJ; 



and OLn = 1, one finds: 



■■,ed t C £ -\^V-A{x)YC £ 
_C £ (z,0) =Co(x) 



V(x)Ce + \Ce\ 2p Ce = in R N x (0, (X)), 



m 



This is the Schrodinger equation with a time-independent external magnetic field. For general facts 
about this equation, we refer again to [Ca] and to the references therein. For the semi-classical 
analysis of standing wave solutions, we refer the reader to the recent work [CJS] and to the various 
references included. For the full (soliton) dynamics, we refer to the recent papers [Se, Sq] which, 
to our knowledge, are the first contributions for this equation. In [Sq], the concentration center is 
precisely the one predicted by the WKB theory. 
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Class III. If m = 1, A = and ctj = jij = Uij = 0, one finds: 

f iedtCe + £ 4&(e ~ V(x)Ce + ^ * \(e\ 2 (e = in R N X (0, Oo), 

lC £ (^0)=Co(x) inR". 

This is the Hartree or Newton-Schrddinger type equation. For basic facts about this equation, we 
refer again to [Ca] and references therein. For the study of standing waves in the semi-classical 
regime, we refer to [WW] and the references included. The physical motivations for these equations 
were detected by Penrose who derived the Schrodinger-Newton equation by coupling the linear 3D 
Schrodinger equation with the Newton law of gravitation, yielding 

iedtCe + ^ AC, - V{x)Ce + ^e*Ce = in R 3 x (0, oo), 
-e 2 A* £ = ^|C e | 2 inM 3 , 
where /i is a positive constant. Of course, this system is equivalent to the nonlocal equation 

IedtCe + " V(x)( £ + * £ * |C £ | 2 C e = in R 3 x (0, oo), * e (x) - 



For the study of point-particle dynamics for this equation with smooth nonlocal potentials, we refer 
the reader to [FTY], where the authors follow an approach different from that used in [BJ, Ke]. 

Class IV. If m = 1 and ctj = jij = Uij = 0, one finds: 

f iedtCe - KfV - A(x)) 2 C e - V{x)Ce + ^ * \Ce\ 2 Ce = in R N x (0, oo), 

\Ce(x,0) = Co(x) inM^. 

This is the Hartree type equation with magnetic field. As for the previous cases, concerning the basic 
facts about this equation, we refer to [Ca]. With respect to the semiclassical analysis of standing 
waves we are not aware of any paper. The soliton dynamics behaviour is contained in the present 
paper for smooth potentials. 

Class V. If m = 2, A = and (3j = u>ij = 0, one finds: 

'iedtd + T A C £ X - V(x)Q + a 1 |Cif P C £ 1 +7i2|C £ 2 r 1 |C £ 1 r 1 = inl w x (0,oo), 
ied 4 C £ 2 + ^ AC| - V{x)<* + a 2 |C|| 2p C £ 2 +7i2|C £ 1 | p+1 |C £ 2 |^ 1 = in R N x (0,oo), 
Ux,0)=( (x) inR N . 

This is the weakly coupled Schrodinger system with two components. With respect to the semiclas- 
sical analysis of standing waves, in the last few years the interest for this systems has considerably 
increased. We refer for instance to [AC, LW, MPS1, Si] for the study of the structure of the associ- 
ated ground states solutions (vector versus scalar ground states depending upon the strength of the 
interaction 712 > 0). For the behaviour in the semiclassical limit, we refer the reader to [DW, MPS1]. 
The soliton dynamics behaviour is contained in [MPS2, MPS3], essentially in the ID case. 

1.2. The main result. In this section we shall provide the suitable background allowing us to 
formulate the statement of the main theorem of the paper. 

1.2.1. Framework and main ingredients. Throughout this paper we denote by Ha,s the Hilbert 
space defined as the closure of C^°(R N ; C m ) under the scalar product 

(u, v)h a e = K / (Du ■ Dv + V (ex)uv)dx , 
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where Du = (D\u, . . . , D^u) and Dj = i~ 1 dj — Aj(ex), with induced norm 



MIL, 



J \vu — A{ex)u dx + J V{ex)\u\ 2 dx < oo. 



The dual space of Ha, £ is denoted by H' A , while the space H\ is the set of u such that 

2 



1 2 



It 



L 1 



+ 



qV- A{ex)) 2 u 



< oo. 



l m ) is equipped with the standard norm ||</>||^i 



IW||£ 2 + II0III2. We study 



(D) 



Finally, 

problem (P) for an initial datum </>q : M. N — ► C m given by 

(/) ^o(x) = rj(a - X£ (o)) e i ^(^(°))-(-^(°))+-^( )], j = 1, . . . ,m 

where xq/e and £0 are the initial position and the initial velocity in HH N of the following first order 
differential system 

'x £ (t)=Ut), 

Ut) = -eVV(ex £ (t)) - e&(t) x B(ex e (t)), 
x £ (0) = f, 

U 6 (o) = &, 

with B = V x i. Notice that, for the solution of (D), we have 

fi(t) = e(t), 

£(t) = -VV(x(t))-Z(t) xB(x(t)), 
x(0) = x , 

U(o)=6>. 

The rescaled components (a;(t),£(i)) of system (1.3) might appear in the proofs of some result. 
Notice that the initial datum referred to the original problem (S) reads as 



(1.3) 



x £ 



jA 



, j = 1,... ,m 



This is the usual formula for the (soliton) initial datum considered in [BJ, Ke] when A = and 
in [Se, Sq] when A 7^ 0. Furthermore, we assume that r = (n, . . . ,r m ) € H l (R N ,R m ) is (up to 
translation) a real ground state solution of the elliptic system 



(S) 



j = 1, ... ,m, 



I |2p 



m 



with respect to the notation of | • L previously introduced. We also set 
(1.4) m» := ||r,-||? a , j = l,...,m, 



Notice that, setting for all t € S 

(1.5) 
where 

A4 



1 



W(0 = ^ e (t)| 2 + ^x £ (t)) + .M, 



£(0)j 
2M 



» m m 
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it follows that TL is a first integral associated with (D), namely 

H(t) = H(0) = ^o| 2 + V(x ) + M, for all t G K+. 
In turn, the function TC is independent of both time and e > 0. 
1.2.2. Assumptions on the potentials. We first give the following 

Definition 1.1. Consider the potentials V : R N ^ R, A : R N — » R N and : R N — » M and 
a ground state solution r of (S) which is chosen to build up the initial datum (I). We say that 
(V,A, 3>,r) is an admissible string for the point particle dynamics of problem (P) if r, is radially 
symmetric, Xirj G L 2 (M Ar ) /or a/Z i = 1, . . . , iV and j = 1, . . . , m and the following Properties 1.2 
(well-posedness) and 1.3 (non- degeneracy /energy convexity inequality) hold true. 

Property 1.2 (Well-posedness). Assume that < p < 2/N . Then, for all e > and 4>o G H^ 6 , 
there exists a unique global solution 

( / }£ eC(R+,H A>£ )nc 1 (R+,H' A>£ ), 

of problem (P) with sup ||</> e (i)||# A < oo. Furthermore, the mass Mi associated with </4(i) , 

Mi{t) := J \4(t)\ 2 dx, t G R + , j = l,...,m, 

and the total energy E £ , 




are conserved in time, namely 

A/|'(t) = A/|(0) and E £ {t) = E e (0), for all t£R + , j = 1, . . . , m. 
Finally if 4> G H\ £ , then <f> £ G C(R+, H\ £ ) n C 1 (R + ,L 2 (R N ; C m )). 

We also consider the functional £ : H 1 (R N ;R m ) — > R associated with system (S) 

£(u) = - ivuix^dx-^^- / \ Uj (x)\^dx- E / K^r'K^r 1 ^- 

In a large range of relevant situations, a ground state solution r of (S) satisfies the characterization 

(1.6) £ (r) = min{£ (n) : n G fT^R^.M" 1 ), ||n|| L 2 = ||r|| i2 }. 

For m = 1 this is a classical fact. For m = 2 see e.g. [MMP1]. 
We consider now the following 
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Property 1.3 (Non-degeneracy/Energy convexity inequality). There exist two positive constants 
C and C such that the following condition holds: if U 6 iJ 1 (R JV ; C m ) is such that ||C/||^2 = ||?"||£2, 
where r is a ground state solution of (S), then 

(1-7) Tu<C(£(U)-£(r)), 

where 

(1-8) T u= inf \\U(-)-(e ie W 1 (-+y),...,e w ™r m (- + y))f Hl , 

y£W. N 
6»i,...,6l m e[0,27r) 

provided that Tu < C 

The energy convexity inequality is essentially a feature of a ground state solution r. It is generally 
a quite delicate issue to consider, based upon nontrivial spectral estimates and the fact that the 
kernel of the linearized operator is ^-dimensional and spanned by the partial derivatives djr of 
r. Let us point out which is the current knowledge of particular cases, within our framework, 
where this assumption is indeed satisfied. For the Schrodinger equation with or without magnetic 
field, Property 1.3 is satisfied, since the (unique) ground state solution of -^Ar + r = r 2p+1 is non- 
degenerate and satisfies suitable spectral estimates (see the striking works of Weinstein [Wei, We2]). 
For systems, already in the case of two components, the situation is still very far from being 
completely understood. On the other hand, very recently Dancer and Wei have proved in [DW] the 
existence of non-degenerate ground state solutions in some particular cases, providing an important 
tool in connection with Property 1.3. In the one dimensional case, Property 1.3 has been verified 
in [MPS3] for two-components weakly coupled nonlinear Schrodinger system. The main obstacle in 
dealing with the higher dimensional case is the smoothness of the energy functional £ which is not 
of class C 2 due to the presence of the coupling terms f |(/>*| p+1 |i^ ? | p+1 , being p < 2/N < 1. 

1.2.3. Statement of the result. On the external potentials V and A, on the nonlocal term $ and 
on the ground state solution r of (5) which is chosen to build the initial datum (I), we assume 
that they are admissible for the point particle dynamics in the sense indicated above and that the 
following conditions hold: 

(V) V £ C 3 (R N ) is positive and \\V\\ C 3 < oo ; 
(A) A € C 3 (R N ;R N ) with \\A\\ c s < oo; 
($) $ € C 3 (R N ) positive with ||$|| C 3 < oo. 

We shall think $ as a smooth function decaying at infinity as \x\~ p for some p > (for instance, in 
M. N with N > 3, decaying as the Coulomb potential |x| 2_Ar ) having a maximum point at the origin. 

Under the previous assumptions, we can state the main result of this paper. 

Theorem 1.4. Assume that $ = in the vectorial case m > 1. Let <j) £ be the family of solutions 
to problem (P) corresponding to the initial datum (I) modelled on a ground state r of (S) and let 
(x £ (t),£, £ (t)) be the solution of (D). Then there exist 5 > 0, Eq > and shift functions 6\, . . . , 0™ : 
K + — ► [0, 2n) such that, if ||-A||c 2 < 8, then 



4(x,t) = e m(tyx+8t(t)+A(ex £ (t)y(x-x £ (t)) r .( x _ Xe ( t )) +w J( a . )t ) > 
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where \\uji(t)\\jji < 0(e), for alls 6 (0,£o) end j = l,...,m, locally uniformly in time with the 
time scale e" 1 . Furthermore, without restrictions on \\A\\ C 2, there exists Eq > such that 

(1.9) \4(x,t)\= rj (x-x £ (t))+u4(x,t), 

where ||a>f < 0{e), for all e E (0, £o) an d 3 = 1, . . . , m, locally uniformly in time with the time 
scale e -1 . 

This kind of results has the origin in some works in linear geometric asymptotics which go back to 
the 70's (see [GS]). We stress that, in the vectorial case m > 1, we are not aware of any physically 
reasonable model including the nonlocal coupling terms. Hence, for m > 1, we consider systems of 
coupled Schrodinger equations with local terms, which are being extensively studied in the literature 
of recent years. 

Remark 1.5. Rescaling back to problem (S), the approximated representation formula reads as 

locally uniformly in time, where we have set l d J e (t) = £#e(^/ e ) an d ^e(x,t) = ui(x/e,t/e), which 
reads as in [Sq] and in the previously cited papers in the particular cases m = 1, A = and $ = 0. 



Plan of the paper. 

In Section 2, we prove various preliminary Lemmas, particularly focused on the asymptotic be- 
haviour of the energy, for e small. In Section 3, we prove some Lemmas, focused on the asymptotic 
behaviour of the density and of the momentum associated with the solution, for e small. In Sec- 
tion 4, we prove a result yielding a precise control on the norm of the error function uP e which appears 
in Theorem 1.4. Finally, in Section 5, we conclude the proof of the main result, Theorem 1.4. 



Notations. 

(1) The imaginary unit is denoted by i. 

(2) The conjugate of any z € C is denoted by z, the real and imaginary parts by and Qz. 

(3) The symbol M + means the positive real line [0, oo). 

(4) The ordinary inner product between two vectors a, b 6 is denoted by a ■ b. 

(5) The standard LP norm, 1 < p < oo of a function u is denoted by ||m||lp. 

(6) The symbols <9 t and dj mean ^ and respectively. A means + • • • + -§^r- 

(7) The symbol C k {R N ; C m ), for k G N, denotes the space of functions with continuous deriva- 
tives up to the order k. Sometimes C k (M N ;C m ) is endowed with the norm 

H0ii c * = Yl iwiu~ <°°- 

\a\<k 

(8) The symbol f f(x)dx stands for the integral of / over M. N with the Lebesgue measure. 

(9) The symbol C 2 * denotes the dual space of C 2 . The norm of a v in C 2 * is 

|M| C 2* = sup{| J 4>{ex)vdx\: <p£C 2 {R N ), \\4>\\ C 2 < lj. 

Clearly, C 2 * contains the space of bounded Radon measures. 

(10) C denotes a generic positive constant, which may vary inside a chain of inequalities. 

(11) 0{e) is a generic function such that the limsup of e~ 1 0{e) is finite, as e — > 0. 
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2. Some preliminary stuff 

Observe that, from Property 1.2, due to the choice of the initial datum (/), the masses Mi (t) 
are also independent of e. Indeed, via the mass conservation law, by the form of the initial datum 
and (1.4), we have 

(2.1) Mi(t) =Mi(0) = J \4(x,0)\ 2 dx = J rj (x-x e (0)j 2 dx = \\ rj \\ 2 L2 = mj , 

for all e > 0, t € R + and j = 1, . . . , m. 

We now recall a useful identity (see e.g. [Ke, Lemma 3.3]). 

Lemma 2.1. Assume that g : R N — > R is a function of class C 2 (R N ), \\g\\c 2 < 00 > ara ^ r is a 
ground state solution of (S). Then, as e goes to zero, for any i = 1, . . . ,m it holds 



g(ex + y)r i (x)dx = j g(y)r i (x)dx + 0(e 
for every y £ R^. 

In a similar fashion, we have the following counterpart to be used for the nonlocal term. 

Lemma 2.2. Assume that g : R N — > R is a function of class C 2 (R N ), ||(?||c2 < oo, and that r is a 
ground state solution of (S). Then, as e goes to zero, for any i, j = 1, . . . , m it holds 



g(e(x — y))r { (x)rj(y)dxdy = mimjg(0) + 0(e ). 
Proof. By Taylor expansion, for some point £ of the form £ = et(x — y) with r 6 (0, 1), we have 
g(e(x - y))r 2 {x)r 2 Ay)dxdy = 



= 9(0) / / r 2 (x)r 2 (y)dxdy + e^D h 5(0) ■ (x h - y h )r 2 {x)r 2 (y)dxdy 
JJ h=1 JJ 

e 2 N ff 

+ y X] // D lk9{0{xh - Vh)(xk ~ y k )r 2 {x)r 2 (y)dxdy 
h,k=l J •* 

N 

= m i m j g(0)+e^2D h g(0) / x h r 2 (x)dx J r 2 {y)dy 

h=l 

-eJ2D h g(0) [ y h r 2 (y)dy f r 2 (x)dx 
h=i J J 

e 2 N ff 

+ Y Yl // D hk9{0{xh - I/a) (a* - y k )r 2 {x)r 2 (y)dxdy 
h,k=l J J 

e 2 N ff 

= mimjg(0) + — ^ / / D 2 hk g{£){x h - y h )(xk - yk)r 2 {x)r 2 {y)dxdy 



h,k=l' 
2\ 



= mimjg(0) + 0(e). 

In the above computations we used the fact that \D\ k g{Cj\ < ll<?llc 2 < °°j that, since r\ is radially 
symmetric, f Zhrf(z)dz = and, finally, that zyji G L 2 (K N ) for any h and i (cf. Definition 1.1). □ 
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In the next result we obtain an asymptotic formula for the energy, linking the functionals E £ , £ and 
TC, up to an error 0(e 2 ) (see also [Sq]). 

Lemma 2.3. For every t G M + , as e goes to zero, it holds 

E £ (t) = £{r) + MH(t) + 0(e 2 ). 
Proof. Taking into account that, in view of Lemma 2.1, for all j = 1, . . . , m we have 
J r 2 (x)\A(ex + x )\ 2 dx = \A(xo)\ 2 mj + 0(e 2 ), 

J r](x)A(ex + x ) ■ (A(x ) + £o)dx = A(x ) ■ (A(x ) + ^) mj + 0(e 2 ), 
as e goes to zero, it is readily checked that, for any j = 1, . . . , m, we get 
J | — (a; — cc e (0)) e i[^(^o) ^-^ e (o)) +;c Co]^ | 2 da - = J \V rj (x)\ 2 dx + mj \^\ 2 + 0(e 2 ). 

In turn, by combining the conservation of energy (see Property 1.2) and the conservation of the 
function TL (see definition (1.5)), taking into account Lemma 2.1 and Lemma 2.2, as e goes to zero, 
we get 

E £ {t) = P e (0) = E £ (r(x - x £ {0))e iWxoHx - x ^ 0))+x < a] ^ 

m „ . _„ 2 



J{A(x )ix-x,(0))+x-Zo] 



= \1L I \{j- A (exj)(rj(x-x e (0))e 

m „ m „ m „ 

+ £ / V(x + ex)r 2 (x)dx-Y,^ l T / |r/ p+2 dx - / Inf+^r^dx 

3=1 3=1 i,j=l : ijij 

"Ey// *M*-v))\rj(*)\ 2 \rM\ 2 tedV-f2 ^ J J ~ vWii^W^dxdy 

3=1 

m ,. m 

= £(r) + Y^ / V(x + ex)r 2 (x)dx + -^2m j \£ \ 2 + MM + 0(e 2 ) 

j=i J 3=1 
m 1 m 

= £{r) + J2 mjV(x ) + g J2 m Mo\ 2 + MM + 0{e 2 ) = £{r) + MH(t) + 0(e 2 ). 

3=1 3=1 

□ 

The function : R N x ]R + — > M. m+N is the (magnetic) momentum of <p £ , defined as 
(2.2) P £ (x,t) := Q(4> £ (x,t)(V(j) £ (x,t) -iA(ex)(f) E (x,t))), x G R N , t G R + . 

Then, we have the following 

Lemma 2.4. Let <f> £ be the solution to problem (P) corresponding to the initial datum (I). Then 
there exists a positive constant C such that 

\\r 1 V4> E {- ) t)-A{ex)4> £ {;t)f L2 <C, 

for all t G M + and any e G (0, 1] . In particular, 



sup | / p e (x,t)dx\ < oo. 



10 



R. SERVADEI AND M. SQUASSINA 



Proof. By Property 1.2 the total energy E e is conserved and, in addition, can be bounded indepen- 
dently of e (due to the choice of initial datum, see Lemma 2.3). Taking into account the positivity 
of V and the definition of E e , it follows that there exists a positive constant C such that 



vV0 e (-,t) " A(ex)<f> e (-,t) 



L- 



-Vcj) e (x,t) - A(ex)<p e (x,t) 



dx 



2 r 

2E £ (t)-2 I V(ex)\<f )£ (x,t)\ 2 dx + — Y,a j / \4(x,t)\ 

p + j= i J 



2p+2 



dx 



+ 



P + 



m „ 



(2.3) 



+ E^7/ He(x-y))\4(x,t)\ 2 \4(y,t)\ 2 dxdy 

3=1 

m „ . 

+ E u l3 JJ <b{e(x-y))W £ {x,t)\ 2 W e {y,t)\ 2 dxdy 

i,3,i¥=3 

3=1 i,3,i¥=3 
m „ „ 

+ E& /y ^ - y)M(*)l 2 l^(y)l 2 ^y 



3=1 
m 



+ E w « J J H<x-y))\4(x,t)\ 2 \4(y,t)\ 2 dxdy. 

i,3,#3 

By combining the diamagnetic inequality (see e.g. [EL] for a proof) 



(j-A(ex) 



a.e. m 



DiV 



with the Gagliardo-Nirenberg inequality, setting $ 



pN 
2p+2 



€ (0, 1), we obtain 



i(;t)\\ L ^ +2 < \\4(-Mi\M4(-,t)\\\l2 < \\ti(;t 



i 

1 1,2 



V 



A(ex))4(;t) 



L 2 



for any j = 1, . . . , m. While, by the conservation of mass, we deduce that 

=Mi(t)= mj j = l,...,m, 
independently of s (see formula (2.1)). Hence, for all e > 0, we get 



(2.4) 



l#(-,t)llg&< 



(i-fl)(p+i) 



1 . pN 

1 



for any j = 1, . . . , m and for some positive constant C, where we have set, for f > 0, 
T e (t) = max T£(t), T>(t) = ||-V^(-,t) - A(ex)^(-,t)|| i2 . 

j=l,...,m 1 

Observe also that, as $ is uniformly bounded, for any i, j = 1, . . . , m we have 



$( e (x-y))|^(x,t)| 2 |^(y,i)| 2 ^dy <C / ^(Mjfdz / |<^(y, t)| 2 (iy = Cm m 
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Finally, notice that, by Young inequality 



(2.5) 



J l^^^)!^ 1 !^^,*)!^ 1 ^ < 1 1|^(-, llf^ + | II 0i II f^+ 2 , 



for any j = l,...,m. Putting now together all the previous inequalities from (2.3) to (2.5), we 
finally obtain (T £ (t)) 2 < C + C(T £ (t)) pN for t > 0. Taking into account that pN < 2 by the 
assumption on p, if T e (t) was unbounded with respect to t or e, the above inequality would yield a 
contradiction. Hence T e is uniformly bounded with respect to t and e, so that the first assertion of 
Lemma 2.4 holds. In order to prove the final assertion observe that, taking into account the mass 
conservation law, by Holder inequality we get 



p e (x, t)dx 



< / \ P ?(x,t)\dx<\\M;t) 



rV0 B (-,t)-A(es)0 e (-,t; 



L 2 



< a 



for all t € M + . The assertion follows by taking the supremum over t in M + . 



□ 



For the next lemma we need to introduce the total magnetic momentum defined as 



q?(x,t) = ^2(p?y(x,t), x € M. N , 

3=1 



Then, on a suitable function tp e (related to the solution <p £ ), we have the following 

Lemma 2.5. Let cp e be the family of solutions to problem (P) corresponding to the initial datum (I). 
Let us set, for any e > 0, t € M + and x £ 



(2.6) 



il4(x,t) = e-'^y^+^h-' 1 ^^ 4 (x + x £ (t),t) , j = l,...,m, 



where (x e (t),£ e (t)) is the solution of system (D). Then, as e goes to zero, 



£{^ £ {t))-£{r) = MH(t)- j V{ex)\ck £ (x,t)\ 2 dx + -M\Ui) + A{ex £ (t))\ 2 



- (£ £ (t) + A(ex £ (t)) ■ j q?(x,t)dx - (&(t) + A{ex £ {t)) ■ J A(ex)\<f> E (x, t)\ 2 dx 
+ \A{ex)\ 2 \4) £ {x,t)\ 2 dx + j A(ex) ■ q£(x,t)dx 

+ \p^j JJ -y)M(x,t)\M(y,t)\ 2 dxd y 

+ \ E ^JjHe(x-y)M(x,t)\ 2 \4( y ,t)\ 2 dxdy + 0(e 2 ). 

Proof. By a change of variable we see that HV'e(i) ||^ 2 = m j f° r j = 1, ... ,7Ti. Hence the mass of 
tp £ (t) is conserved through the motion. Let pi(x,t) = ^s((j) J £ (x, t)V(fr{(x, t)j for x G R N , t € M+ and 
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j = 1, . . . , m be the j-th magnetic- free momentum. A direct computation yields 

If I m 

£mt)) = - / \VMx,t)\ 2 dx + -Y,m j \Ut) + A(ex e (t))\ 2 

J 3=1 
m „ -. m „ 

-^(Ut) + A(ex £ (t)). p>(x,t)dx-—J2<X3 / We{^t)\ W d2 
j=l J P+ j=l J 

m 

— E Hi] \4(x,t)\ P+1 \4M\ P+1 dx, 



p+ - 

h3, *7=? 



so that we obtain 

1 r 1 2 

f/.r 



1 2 

T V0 e (x,i) - A(ex)^ e (z,t) 
i 



J \A{ex)\ 2 \<j) £ (x,t)\ 2 dx + f2 J A(ex)-pi(x,t)dx 
+ jE^'I^C*) + ^ £ (t))| 2 - E(£ B (t) + A(ex 6 (t)) • / pl(x,t)dx 

m*,t)\ 2p+2 dx-— -raj m*M p+1 \4M\ v+1 dx. 

i=i 



Then, taking into account the definition of and of 7i and Lemma 2.3, we obtain 

£{^ £ (t)) - 8(r) = MH(t) - j V{ex)\<t> £ {x,t)\ 2 dx + ^M\£ e (t) + A(ex e (t))\ 2 



/m 
Y,viMdx 

+ ^ E ^JJ$(e(x-y)M( X ,t)\ 2 \rt(y,t)\ 2 dxdy + 0(e 2 ), 



i,3, i^j 

as e goes to zero. Finally, since pi(x,t) = (pfy(x,t) + A(ex)\(f>i(x, t)\ 2 and recalling the definition 
of qf, we obtain the desired conclusion. □ 

Now let us introduce two functionals in the dual space of C 2 

(2.7) j U l £ (x,t) ■ cp(x)dx = j ip(ex) ■ q*(x,t)dx - M<p(ex e (t)) • £ e (t), V<^ e C 2 (R N ; R N ), 

(2.8) [ U 2 (x,t)ip( X )dx = [ ip(ex)\^ £ (x,t)\ 2 dx- Mip(ex e (t)), V92 € C 2fiaN ■ m ^ 
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for all t £ R + , where M is given in formula (1.4). Moreover, define the function O e : 
n E (t) = n E (t) + pf(t), where 

(2.9) f2 e (t) := y n^(a;,t)(ix 

and 



13 
as 



+ sup 




+ he(t)\, 


t G M + 


HV|| C 3<1 









t € 



Ul(x,t) ■ A(x)dx , 
%(t) := Mex £ (t) - / sxx{ex)\(j) £ {x,t)\ 2 dx, 



t e M + , 

where % € C7°°(R Ar ) is such that < x < 1, x( x ) = 1 in B(0,p) and xi x ) = in 
a suitable /5 > that will be suitably chosen later. 

Now we are able to prove an estimate on the energy of ift e . 



\B(0,2p), for 



Lemma 2.6. Assume that $ = if m > 1 and let ip £ be the function defined in formula (2.6). 
Then there exists a positive constant C independent of e such that 

< £(Mt)) - £(r) < CSl e (t) + 0(e 2 ), 

for all t £ M + and any e > 0. 

Proof. We claim that S1 £ (0) = 0{e 2 ) as e goes to zero. In fact, by definition of Q £ , we have 



(2.10) 


Q e (0) = 




+ sup 










IMI C 3<1 





First of all, let us estimate the first term in the right-hand side of (2.10). Taking <p = 1 in (2.7) 
and using (I), we get 

n;(s,o)^ = yV(*,o )<fc -M f (o) 

= £ / ^(^(x.OJCV^Cx.Oj-LlCe^^Cx.Oj^dx-Afeo 
7=1 7 

m 



7=1 



/ r?(x - x e (0)) A(x ) + Co - A(ex)] cfo - M£ 



M^4(x ) - / r i( x ~ Ze(0))A(ea;)cfa 

7=1*' 

MA(x ) - Y / rf(a:)A(ex + xo)da; = C(e 2 



as e goes to zero, in light of Lemma 2.1. In a similar fashion, one gets pf(0) = 0(e 2 ). Now consider 
the second term in the right-hand side of (2.10). Let ip E C 3 (R ) with |Jy ||c? 3 — !■ Then, 



n2(.,oM^ = /«. (l .o)|^-M, W o>) 

m „ 

= ^2 vfa + x )r](x)dx - Mip(x ) = 0{e 2 ) 

.7=1 ^ 
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as e goes to zero, again using Lemma 2.1. We finally estimate 7 e (0). As above we have 



7 e (0) = Mx(0) - e J x X (ex)\(f> £ {x,0)\ 2 dx 

TTL n III/ j\ 

= Mxq — / xx(ex)r 2 (x — x e (0))dx = Mxq — / (ex + xq)x(sx + xo)r 2 (x)da 
= Mx -^ / x x(x )rj(x)dx + O(e' 



Ms (1 - x(a?o)) + 0(0, 



thanks to Lemma 2.1. Now, from [Kc, Lemma 3.1-3.2] (where one has to use the 5 a at some point a 
is defined as (<5 , (p) = (p(sa) for all <p 6 C 2 (M N )), we learn that there exist three positive constants 
Kq, K\, K2 such that, for all y,z G E , K\\ey — ez\ < — <5 Z || C 2» < K<i\ey — ez\, provided that 
\\5 y — 8z\\c 2 * ^ Ko- Let then p = K\ sup e6 r u sup te r 0j y / E i |ea: e (t)| + i^Oj where To > is fixed (to be 
chosen later on, see Lemma 3.4). Then, in view of the definition of x, we obtain that 7 e (0) = 0(e 2 ) 
as e goes to zero, since \xo\ < p. Hence the claim is proved. 

Now we are ready to prove the assertion of Lemma 2.6. By using Lemma 2.5, the definition of 
H, (2.7) and (2.8) we obtain 



£(Mt))-£(r) 



1 



M%[t)\ 2 + MV{ex £ {t)) + MM 



1 



V(ex)\(f) £ (x,t)\ 2 dx + -M\Ut) + A(ex £ (t))\ 2 



Ul(x,t) (Ce(t) + A(ex e (t))) 



dx-M 



Ut) + A(ex £ (t)) • 



- (f B (t) + A(ex £ (t)) ■ ( J n 2 (x, *)A(a;)ds + MA(ex e (t)) N 
+ |A(ex)| 2 |0 £ (x,t)| 2 dx + y n^(x,t)yl(x)(ix + MA(ex e (t)) -^(t) 

+ fh{e{x-y))W e (x,t)\ 2 \^{ y) t)\ 2 dxdy 



i=i 

- y 
2 ^ 



Uij a $(e(x-y))\4(x : t)\ 2 \4(y,t)\ 2 dxdy + 0(e 2 ). 



Let us set (with the convention that ton = Pi 



LOij 



f <S>(e(x-y))\<t>l(x,t)\ z \(tP £ (y,t)\ 2 dxdy-<S>(0) ^ ^m l m ] 

i,j=l i,j = l 
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In turn, using the definition of M, we have 
£(Mt))-£(r) 

< n e (t) + MV(ex £ (t)) - J U 2 e (x,t)V(x)dx - MV{ex e {t)) + ^M\A(ex £ (t))\ 2 

-J l£(x,t)[&(t) + A(ex&))]dx - &(t) + A(ex e (t)) f U 2 £ (x, t)A(x)dx - M\A(ex £ (t))\ 2 

+ ^Ju 2 (x,t)\A(x)\ 2 dx + ^M\A(ex £ (t))\ 2 + J ' U l £ (x,t)A(x)dx + 0(e 2 ) 

= Ve(t)~ J U 2 (x,t)V(x)dx- J U l £ (x,t) [&(i) + A{ex £ {t))]dx 

- (&(t) +A{ex £ {t)) j n 2 £ (x,t)A(x)dx 

+ U 2 £ (x,t)\A(x)\ 2 dx + j U 1 £ (x,t)A(x)dx + 0(e 2 )< 7]E (t) + Cn £ (t) + 0(e 2 ), 

for e sufficiently small. If m > 1 we assume that <I> = 0, and the assertion follows. If instead m = 1, 
observe first that from definition (2.8), by choosing <p(x) = &(x — ey) and tp(y) = $>(ex e (t) — y) 
respectively, we have 

j &(ex - ey)\4>\(x,t)\ 2 dx = j Tl 2 (x,t)&(x - ey)dx + mi$(ex £ (t) - ey), 
mi J <$>{ex £ {t) - ey))\4> l £ (y,t)\ 2 dy = m x J U 2 (y, t)$(ex e (t) - y)dy + *(0)m?. 



In turn, we have 
Ve 

= C 



(t)<C J J n 2 (x,t)<$>(x - ey)dx + m 1 <$>(ex £ (t) - ey) 



\<l>l(y,t)\ 2 dy - *(0)m? 



n 2 (x,t)$(x - ey)dx 



< Crri\ sup 

IMIcS^ 1 



< C sup 

IMI C 3<1 



H 2 (x, t)(f(x)dx 



l(y,t)\ 2 dy + mi / *(ex e (t) - ey)\4> l £ {y, t)\ 2 dy - $(0)m? 



+ Cmi 



U 2 (y,t)^(ex £ (t)-y)dy 



t)(p(x)dx 
In turn, we conclude that 



< cn £ (t) < cn £ (t). 



£{^ £ (t))-8(r)<CSl £ {t) + 0{e 2 ) 
as e goes to zero, for some positive constant C. Hence the proof of Lemma 2.6 is complete. □ 

Since the function {t ^ $7 e (i)} given in (2.9) is continuous and recalling that fi e (0) = 0(e 2 ) as 
£ — > (see the proof of Lemma 2.6), for any fixed To > and uq > 0, we can define the time 

(2.11) T* := sup{t G [0,T o /e] : Q e ( a ), I> e(s) < <t , for all s e (0,t)} > 0, 

for any e > 0, where is defined according to (1.8) and r^ e ( ) = 0. Now we are able to provide the 
main result of this section, related to a representation formula for the solution <fi £ of problem (P). 
For the proof, it is enough to adapt the proof of [Sq, Theorem 4.2]. The fact 

Theorem 2.7. Let (f> £ be the family of solutions to problem (P) corresponding to the initial da- 
tum (I) modelled on a ground state solution r of problem (S) and let (x £ (t),( i£ (t)) be the global 
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solution of (D) . Then there exist positive constants eo and C , locally bounded functions 8j, . . . , 6™ : 
R + -> [0, 2vr) and y £ : R + -> R N such that 

4(x,t) = e KUt>x+em+A(ex s (t)>(x-*At)) r .( x _ y e ( t ))+<4(t), 



where \\u J £ (t)\\ H i < Cy/Q £ {t) + 0(e), for all e G (0,e ), t € [0,T*) and j = 1, ...,m. 

3. Density and momentum identities 

This section is devoted to some important identities involving the momentum p^ and the total 
magnetic momentum qf related to problem (P). 

Proposition 3.1. Let (j) £ be the solution to problem (P) corresponding to the initial datum (I). 
Then the following identities hold true 

d\4\ 2 

(3.1) -^(x,t) = -div x (pAy(x,t), x£R N ,t£R + ,j = l,...,m, 

/ ~^^ x ^ dx = ~ J 1^(x,t) X eB(ex)dx - J eW(ex)\(f) £ (x, t)\ 2 dx 

(3 . 2) JJ £VHE(x-y))\4(y)\ 2 \4(x)\ 2 dxdy 

3=1 



u i:j J j eV$(e(x - y))\4(y)\ 2 \4(x)\ 2 dxdy, 



i,3=l,i¥=j 

for t € R + , where B = V x A is the magnetic field associated with A. 

Proof. The proof follows the lines of the corresponding proof in [Sq] for the scalar case without the 
presence of nonlocal potentials. By formula (2.2), for any j = l,...,m, (p^y is the vector whose 
components, which we denote by (p^)}, are given by (p^)l = Q{^s(x, t)(d£(pi(x, t)—\Ai{ex)(jy £ {x, i))) , 
for £ = 1, . . . N. Let us fix j = 1, . . . , m. Hence 

N 



-div x (pAy( x ,t) = -Y,Hd £ ft E {x,t)(d e 4(x,t)-iA £ (ex)4{x,t)) 



N 



- ^{4M t)(dl4(x, t) - id e Ai(ex)4(x, t) - iA £ (£x)d e 4(x, t)) 

= 2A{ex) -K(V^(x,t)(^(x,t)) - %{4{x, t)A<g(x, t)) + div x A(ex)\4(x, t)\ 2 . 
Moreover, using (P) and taking into account the definition of La, we get 



(x, t) = 23 (fe(x, t) (L A 4(x, t) + V(ex)4(x, t) - \</> e (x, *)|f *) 
-^(ex)*\4> £ \)4{x,t))) 
= -Q(4(x,t)A4(x,t)) + 2A(ex) ■ M(4(x, t)V#(x, t)) + div x A(ex)\4(x, t)\ 2 , 
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so that identity (3.1) holds true. Now let us prove the second one. By definition of the total 
magnetic momentum for any I = 1, . . . , N, we have 

^§r = E ^fr = E (W^) + ^(^))) - E » - E ^ 

3=1 3=1 3=1 3=1 



d\<P: 



mm m 
3=1 3=1 3=1 



and so, integrating over K , it is easy to see that 



(3.3) J ~^dx = Q{dt^ e d t 4)dx " E / 9 ( d e(fedt4)) dx " E / 



3=1 3=1 3=1 



(9/ 



Let us consider the first term in the right-hand side of (3.3). Conjugating the equation, multiplying 
it by 2id^<j) e , £ = 1, ■ ■ ■ ,N, and taking the imaginary part, we have 

2^{d t ^Jd l 4) = -R(A$Jd e 4) + 2A(ex) ■ 3(V0&#) + \A(ex)\ 2 ^ ! d e 4) 
+ div x A(£x)S(4de4) + 2V(sx)U(^J 



2»(|0 B |f $Jd t 4) - 2K ((<D(ex) * \0 £ \ 2 )4d e c%) 



i? 



i=i j=i \ ~ / 

+ 2A(ex) ■ 3?(V<^c^) + \A(sx)\ 2 ^Jd £ 4) 

+ dw x A(ex)^ £ d^ £ ) + d e (V(ex)\<Pl\ 2 ) - ed e V(ex)\4\ 2 

m 



-— j -d i {m 2p+2 )-2 ]r Tiii^rv.r'^v^) 

P i=lMj 

m 

-2p j to((${ex)*\4\ 2 )4id e 4)-2 ]T c^-St ((*(ex) * \4\ 2 )4dt4) ■ 

Hence, integrating over M. N and using the i? 2 -regularity of the functions involved, for all I = 1, . . . , N 
we obtain the following identity 



2 / Q(d t <t> E de4)dx = 2 J A{ex) ■ $X(V<P E dt<P e )dx + J \A(ex)\^((f) £ J d e 4)dx 

+ J dw x A(ex)Q(4dt4)dx - e J d(V(ex)\4 e \ 2 dx 

m „ 

-2 ^ 7ij j mp+^p-^jo^dx 



7|2\ 



2/3, y ft((*( eiE )*|#| 2 )^>)dz 

m „ 

2 £ Wy / R ((*(ex) * \4\ 2 )4A4e) dx. 

i=l,*^3 
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Notice that 



1 

v + 



div x A(ex)Q(ft £ d e 4)dx + 2 A(ex) • 3?(V0£6^)cte = ^E / 9 * Ai 

^ i=i ^ 

Moreover, thanks to the regularity of 4>i, we have 

f; [m^ i \4r i ^jd i 4)dx= f; j i^rv.r 1 ^ ( 4^ ) dx 

i m r 

+1 E y i^r^d^r 1 )^ 

1 m r 

+i E y ^(i^r i i^r i )^ = o, 



i,j = l,i<j ' 



and 



3? ((*(ex) * |#| 2 )0&#) = yy $(e(x - y))|^(y)| 2 K(^(x)a,^(x))^x 

$( £ (x-y)M(y)| 2 c>, [ MMJ! ) dydx 



2 

= ~ JJedMe(x-y)M(y)\ 2 \4(x)\ 2 dxdy 
for all ^ = 1, . . . , N. While, with the same arguments, we get 

u ij Jft{($(ex)*\4\ 2 )4d i 4)dx = ~ f; Uij JJ edMe(x-y)M(y)\ 2 \4(x)\ 2 dxdy. 
Hence, it is easy to see that 



2 J Z(d t e j d e 4)dx = J2f ed e Mex)S{PA4)dx + J \M^)\ 2 d e f-^-J dx 

^y(ex)|^| 2 dx - # JJ ed t *(e(x - y))\4(y)\ 2 \4(x)\ 2 dxdy 
m „ „ 

2 u; tJ II edM^-y)M(y)\ 2 \4(x)\ 2 dxdy 

m „ m „ 

/ e^ i (ex)9(^0i)dx + ^ / eA(^)^Ai(^)|^| 2 dx 



(3.4) i=1 '^ 



i=l i=l 

^y(ex)|^| 2 dx " /3j // ea^(e(s - y))|^(y)| 2 |^(x)| 2 dx f iy 
" E ^ JJ edMe(x-y)M(y)\ 2 \4^)\ 2 dxdy 

for all £ = 1, . . . , iV. As for the second term in (3.3), using again the regularity of (pi, we get 
j Q (d^{^P £ dt4e) \ dx = 0, for any £ = 1, . . . , N. Finally, as for the third term in the right-hand side 



POINT PARTICLE DYNAMICS FOR GENERAL NLS PROBLEMS 19 

of (3.3), by (3.1) we get 

d\4 



h{\ 2 



A t (ex) "^ (x,t)dx = - I A e (ex)dw x (p^) j (x,t)dx 



= J2 / ediMexXp^iix,^ 
i=i ^ 

m „ 

(3.5) = ^2 / zdiA e (ex)Q(ft e (x,t) (difiiix^) - iAi(ex)4>i(x,t))) 

m 

ec?j^(ex)3 : (^(x, t)di4> J e (x, t))dx 

i=i 

m 

y~] e J 4j(ex)9j^(ex)|^(x, t)| 2 <£c 



i=l 

for any £ = 1, . . . , iV. Then (3.3)-(3.5) yield 

a(^) 



at 



(x,t)dx = / e (c^-A^ea;) — c^y^ea;)) ^s{(fP e di^)dx 

m „ 

+ E / £ M £X ) (deAi(ex) - diA e (ex)) \4\ 2 dx 

i,3=l 

m ,. m ,. „ 

/ ea^(ex)|^| 2 dx-J3^ ed e $(e(x-y))\4(y)\ 2 \4(x)\ 2 dxdy 

3=1 J 3=1 J J 

^ wy yy e ^$( e (x-y))i^( y )iv £ (^)i 2 ^ 



q^(x,t) x eB(ex)) £ dx — I edgV {ex)\4>A 2 dx 



- E & yy - y))i^(y)i 2 i^(x)i 2 ^ 

" E "ij JJ edt^e(x-y)M(y)\ 2 \4(x)\ 2 dxdy 

i=l,i^3 

for any £ = 1, . . . , N, so that (3.2) is proved. □ 
Remark 3.2. Taking into account the definition of qf, by (3.1) easily follows 

£ ' (x,t) = -dxv x qA(x,t), x£R N ,t£R + , 



at 

which is consistent with the conservation's laws for the nonlinear Schrodinger equation. 

We now give some estimates on the momentum p^ and the total magnetic momentum qf related 
to problem (P). 

Lemma 3.3. Let <j) £ be the solution of problem (P) corresponding to the initial datum (I) and let 
(x e (t),£ e (t)) be the global solution to (D). Then, in the notational framework of Theorem 2.7, there 
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exist eq > and C > such that 

\\\4P £ (x,t)\ 2 dx - m 3 5 ye{t) \\ c2 , + \\qf £X ^(x,t)dx - M£ e (% e (t) || C 2* < CSl e (t) + 0(e 2 ), 
for every t G [0, T*) and e G (0, Sq) and for all j = 1, . . . , m, where T* is given in (2.11). 

Proof. For any v G H 1 (R N ), we have the formula |V|t>|| 2 = \\7v\ 2 — ^j?^ . Then, by virtue of 
Lemma 2.6, it follows that 



o < em) - £(r) + \Y,j ** < cn E (t) + o(e 2 ) 



5=1' 



for all t G ]R + and e > 0. Moreover, since HIV'sHIl 2 = Ikjlli 2 f° r all j = 1, . . . , m and £(|^ E |) > £ (r) 
by means of (1.6), we have 

(3.6) y -j- < Cn e (t) + 0(e ), 

for every t G IR + and e > and for all j = l,...,m. Following the blueprint of [Sq, Lemma 6.1], 
we get the assertion (see also [MPS2]). □ 

Lemma 3.3 allows us to prove the following result on the distance between the point y £ {t) found out 
in Theorem 2.7 and the trajectory x e {t). For the proof, follow the blueprint of [Sq, Lemma 6.3]. 

Lemma 3.4. In the notational framework of Theorem 2.7 there exist £q > and Tq > (cf. the 
definition ofT* = T*(Tq)) such that 

Px.® - <W)Hc 2 * ^ C\Ex £ (t) ~ ey £ {t)\ < Cn £ (t) + 0(e 2 ), 
for all £ G (0,e ), t G [0,T*), where T* defined as in (2.11). 

Next, we state a strengthened version of Lemma 3.3, obtained thanks to Lemma 3.4. Follow the 
blueprint of [Sq, Lemma 6.4] for a proof. 

Lemma 3.5. Let Tq be as in Lemma 3.4- Let (f> £ be the family of solutions to problem (P) with 
initial datum (I) and let (x £ (t),£ e (t)) be the global solution of (D). Then there exist eq > and 
C > such that 

\\\4(x,t)\ 2 dx - mj5 Xe{t) \\ c2t + \\qf{x,t)dx - M^ £ (t)5 Xe{t) \\ c2t < Cfi 6 (i) + 0(e 2 ), 
for allee (0,e ), t G [0,T*). 

In particular, by the definition offl £ , there exists 5 > with 

(3.7) || \4(x, t)\ 2 dx - mjS^ 11^. + \\qHx, t)dx - M£S¥x s {t) || C 2» < CQ e (t) + 0(e 2 ), 
for all e G (0, £o) an d t € [0,T*) ; provided that ||^4||c 2 < $■ 

Remark 3.6. In Lemma 3.5, while the C 2 *-norm control holds on U J £ = \4> £ {x, t)\ 2 dx — mj5 Xs ^ for 
each j = 1, . . . , m, the control on the momentum holds for the total magnetic momentum qf{x, t). 
This is in fact natural, since looking at the second identity in Proposition 3.1, it is clear that it 
cannot hold for each individual (p^y , unless some other (disturbing) integral terms are added to 
the formula. 
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4. Uniform estimation of Q £ 

Before proving the main result we give an estimate showing that the quantity f2 e (t) can be made 
small at the order C(e 2 ), uniformly on finite time intervals, as e goes to zero. 

Lemma 4.1. Let Tq be as in Lemma 3.4 and £q,5 as in Lemma 3.5. Then there exists C > such 
that h £ (t) < Ce 2 , for all e G (0,e ) and t G [0,T*). 

In addition, if we assume that ||A||^2 < 5 for 5 > sufficiently enough, then Q E (t) < Ce 2 , for all 
e G (0,e ) and t G [0,T*). 

Proof. By the definition of ILi, Lemma 3.5, Proposition 3.1 and system (.D), we obtain 



^-Ul(x,t)dx 
dt eK ' J 



^(x,t)dx - MUt) 



\x,t) x eB(ex)dx + J eW(ex)\cp £ (x, t)\ 2 dx 

III: f r 

+ £ ^« x - y)W e (y)\ 2 We(x)\ 2 dxdy 



3=1 



m „ „ 

+ E u l3 JJ eV§(e{x-y))W £ {y)\ 2 We(x)\ 2 dxdy 

- AfeVV(ea: e (t)) - Me^(t) x B{ex £ {t)) 

If m > 1, we do not have to manage the nonlocal terms, since <3? = 0. If instead m = 1, recalling 
that V$(0) = 0, by Lemma 3.5 and arguing as at the end of the proof of Lemma 2.6, we get 



(4.1) 



eV$(e(x - y^Hy^Kx^dxdy < e[CCl e (t) + 0{e 2 )] , 



for some positive constant C, for all e G (0,£o) an d t G [0, T*). In turn, it holds 
d 



(4.2) 



dt 



U](x,t)dx 



< £ 



+ e 



qf{x,t) x eB(ex)dx + J eW(ex)\(f> £ (x, t)\ 2 dx 

J MeVV(ex)5 Xe{t) dx - j Me£ £ {t) x B(sx)S x ^ t) dx + e[CCl £ (t) + 0{e 2 )} 

j (<£{x,t) - M£ £ (t)8 Xe{t) ) x B(ex)dx 

J VV(ex)(\^ £ (x,t)\ 2 - M5 Xe (t) ) dx | + e [CQ £ (t ) + 0(e 2 )] 
< Ce||^(x, i)dx - Me e (t)^ e(t) || c2 . + Cs\\ \<f> e (x, t)\ 2 dx - M5 Xe[t) \\ c2 * 
+ e [Ch £ (t) + 0(£ 2 )] < e [Cn £ (t) + 0(£ 2 )] , 
for all e G (0,£q) and t G [0, T*). Hence, recalling that O e (0) = 0{e 2 ) as e goes to zero, 



(4.3) 



H\{x,t)dx 



< 



U l £ (x,0)dx 



+ 



I'M ^ l(x ' r)ix 



dr < Ce\l + et) + Ce / S\(r)dr 



22 



R. SERVADEI AND M. SQUASSINA 



for all e G (0,e ) and t G [0,T*). Now, let <p G C7 3 (R JV ) such that \\tp\\ C 3 < 1. Again in light of 
Proposition 3.1 and Lemma 3.5, we have 

(4.4) | J jU 2 £ (x,t)ip(x)dx\ = | J ip(ex)-^- t \Mx,t)\ 2 dx - MeVip(ex £ (t)) ■ &(*) 

- J ip(ex) div x qf(x,t)dx - MeVip(ex £ (t)) ■ 

J eV(p(ex) ■ qf{x,t)dx — j MeVip(ex) ■ ^ £ (t)5 Xe ^dx 

J eVip(sx) ■ (<£(x,t) - MZ £ (t)5 Xe(t) )dx 
< Ce\\qA( x ,t)dx - Mt £ (t)6 Xe{t) \\ c2 » < e[cCl e (t) + 0{e 2 )] 
for all e G (0, £q) and t G [0, T*). Thus, arguing as above, we get 



(4.5) 



sup 

IMI 3<1 



U 2 (x,t)(p(x)dx <Ce 2 (l + et) + Ce / n £ (r)dT 



for all e G (0, £q) and t G [0, T*). Finally, again via Proposition 3.1 and Lemma 3.5, we have 



(4.6) 



\%(t)\ = Me£, £ (t) + J exx{ex) div x q £ (x, t)dx 
= Me£ e (t)- / X7(exx(ex)) ■ q £ (x,t)dx 



V(xx(ex))M£ e (t)8 Xe{ t)dx - J V(xx(ex)) ■ q^{x,t)dx 
< eC\\qA( x ,t)dx - M£ £ (t)5 Mt) \\ c2t < e[cCl £ (t) + 0(e 2 )], 

\ l£ (t)\ < Ce 2 (l + et) + Ce f Cl £ (r)dr, 



which implies 
(4.7) 



for all e G (0, £o) and t G [0,T*). Collecting the above inequalities, recalling the definition of Cl £ (t) 
and taking into account that, for t <T*, by the definition of T* it holds et < eT* < Tq, we get 



Q £ {t) < Ce z (l + et) + Ce / £ (r)dr < Ce* + Ce / n £ {r)dr 

Jo Jo 

for all £ G (0, £o) and t G [0,T*). Hence, Gronwall Lemma yields 

n e (t) < Ce 2 e £t < Ce 2 , 

for all £ G (0, £o) and t G [0, T*), which gives the assertion. Finally, concerning the last assertion of 
the Lemma, recalling again Lemma 3.5 and taking into account the definition of pf(t), if ||^4||c 2 < ^ 
for 5 > small enough, we conclude the proof. □ 
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5. Proof of Theorem 1.4 completed 



5.1. First conclusion of Theorem 1.4. Let To be as in Lemma 3.4 and eq,S as in Lemma 3.5. 
By Lemma 4.1 and the definition (2.11) it follows that T* = Tq/e, for all e G (0, Eq). Hence, 
Q e (t) < Ce 2 for all e G (0,e ) and t G [0,T o /e], in light of Lemma 4.1. Mor eover, by Theorem 2.7 

+ -> [0, 2vr) and y £ : M+ -> such that 



there exist functions 9~,... } 6™ 



where \\u>i 



< Cy/Q e (t) + 0(e), and hence, we have ||(4(t)||Ei e ^ c> ( e )i for a11 e G (0,e ), 



t G [0, To/e] and j = 1, . . . ,m. Lemma 3.4 and Lemma 4.1 also yield \x £ (t) — y £ (t)\ < 0(e), for all 
e G (0,e ) and t G [0,T /e]. Finally, using (A), (V) and (1.5), we get 



(5.1) 



M^t)-x+ei(t)+A(e Xe{ t))ix-x e (t)) (r . (x _ ye{t)) _ r . {x _ Xe(t))) 



< 
+ 
+ 



&(t) + A{Ex e {t)) rj (x - y e (t)) - Tj{x - x&)) 



dx 



Vrj{x-y&))-VTj{x-x e {t)) 



dx 



Tj{x - y £ {t)) - rj {x - x e {t)) dx < C\x e {t) - y e {t)\ 2 < CO{e 2 ), 



for all e G (0,e ), t G [0,T /e]. Therefore, it follows that 



(5.2) 



[( X)t ) - e ^t)-*+ei(t)+A(ex e (t)Hx-x e (t)) r .( x _ Se(t) j 2 ^ < 0(£ 2 ); 



for all e G (0,eo), t G [0, Tq/e] and j = 1, . . . ,m. Hence, Theorem 1.4 holds true in [0, Tq/e]. Now, 
let us take x\ = x £ (Tq/e) and £i = ^ e (To/e) as new initial datum in system (D) and the functions 

4(x) = Tj{x - x eylA(exiy{x-xl)+x^}^ x G M. N , j = 1, ... ,171 



as new initial datum for problem (P). Arguing as above, we can show that Theorem 1.4 holds true 
in [To/e,2To/e] and so, in any finite time interval [0,T/e], with T > 0. The proof of Theorem 1.4 
is now complete under the assumption that ||-A||c7 2 < <5- 



5.2. Second conclusion of Theorem 1.4. To prove the second part of Theorem 1.4, namely 
formula (1.9), we follow the argument of [Se] (which is based upon the original paper by Bronski- 
Jerrard [BJ]). Let us give a brief sketch of the proof. Based upon the identity (see for instance [Se, 
p. 2571]) holding for all v G H l ( 



X7v 



A(ex)v 



\p A ( ex \v)f 



+ IVM 



p A (v) ■= %(v(y v (x, t) - iA(ex)v(x, t))) . 



the energy functional of the Schrodinger problem is rewritten as 

E £ {t) = E?°\t) + E*(t) + E k e {t) + Ef(t), 
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where we have set 

£P ot (i) 



V(ex)\<p £ (x,t)\ 2 dx, 



1 m „ ^ m 

E h M :=-£ / |V|0i|(x,t)| 2 --— Y, 

3=1 J 1 3=1 



a ; 



\4(x,t)\ 2p+2 dx 



m . 

E mx,t)\ p+1 \4M\ p+1 dx, 



ij,i¥=3 

'■ \( P A ( £x \x,t)y\ 2 



\</P e (x,t)f 



3=1 J 

E ?(t) : =-^E^'// He(x-y)M(x,t)\ 2 \4(y,t)\ 2 dxdy 



2 E w « //*(^-y))l^.t)[ 2 l^(i/,t)l 2 ^ 



Notice that, with respect to our notations, we have E^(r\, . . . , r m ) = £(n, . . . , r m ) since are real 
valued and positive functions. Moreover . . . , |if |) = ^(|^|, ■ ■ ■ , \(f>?\) = €(\4>l\, ■ ■ ■ , \4>f\)- 

At this stage, keeping in mind that we possess Lemma 2.3, which expands the energy E e (t) up to 
an error 0(e 2 ), by repeating the steps of the proof of [Se, Lemma 3.5], it is readily seen that, as e 
goes to zero, 

< E*(\£\, . . . , |<0 - E h £ (r l} . . . ,r m ) < cCl E (t) + 0(e 2 ). 

This conclusion plays the role of Lemma 2.6 and, as a consequence, by the non-degeneracy/energy 
convexity property (applied with U = (\(j>l\, ■ ■ ■ , see e -§- PJ> Proposition 1] for the scalar 

case), yields 



(5.3) 



e I ) ■ ■ ■ ) I Ye 



M- + ife(t)), • • • , r m (- + y e {t))) Whi < Ctt £ (t) + 0(e 2 ), 



for some y s (t) € 

Moreover, again by the steps of the proof of [Se, Lemma 3.5], we get 



(5.4) 



o<j£(t)--£ 



1 m \f(p A ( £X )( X ,t)y 



3=1 



rrin 



< CCl e (t) + 0(e 2 ), 



as e goes to zero. To achieve this conclusion, one also needs to take into account the following 
elementary inequality (following from the standard Cauchy-Schwarz inequality) 



r a 2 ™ 

/ qf(x,t)dx <M^ 

J 3=1 



f(p A ( ex \x,t)ydx 



t G 



m : 



Furthermore, for any j = 1, ... ,m we have the inequality (see [Se, inequality below formula (28)]; 
see also [BJ, formula (3.2)]) 



{ P A( - £x \x,t)y {J(p A{£x) (x,t)y] 



\4(x)\ 



m*)\ 



dx < 



i r | ( P A ^ (x, t)y\ 2 i | j(p A( -^ ( x , t)y | 



-dx~ 

\4(x)\ 2 2 
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Summing over j = 1, . . . , m, we get 

( P A( - ex \x,t)y (j(p A{£x) (x,t)y) 



2 k 



We{x)\ 



2 1 m 

\4(x)\ dx<E^{t)--Y, 



f{p A( . sx \x,t)y\ 



3=1 



rrin 



In turn, in light of (5.4), we obtain 
(5.5) 



\4i(x)\ 



m, 



\4{x)\ dx < Ch £ (t) + 0(e 2 



as e goes to zero, for any j = 1, . .. , m. Inequalities (5.3) and (5.5) are precisely what is needed in 
order to prove (3.7) of Lemma 3.5 (see the proof of Lemma 6.1 in [Sq], in particular formula (6.5) 
therein; see also the proof of Lemma 4.3 in [MPS2]). Once inequality (3.7) of Lemma 3.5 holds true 
the rest of the proof continues as before, yielding the assertion from inequality (5.3). 
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